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While it is well-known that the electron-electron (ee) interaction cannot affect the resistivity ol a 
Galilean-invariant Fermi liquid (FL), the reverse statement is not necessarily true: the resistivity of a 
non-Galilean-invariant FL does not necessarily follow a behavior. The behavior is guaranteed 
only if Umklapp processes are allowed; however, if the Fermi surface (FS) is small or the electron- 
electron interaction is of a very long range, Umklapps are suppressed. In this case, a term 
can result only from a combined-but distinct from quantum-interference corrections- effect of the 
electron-impurity and ee interactions. Whether the term is present depends on 1) dimensionality 
[two dimensions (2D) vs three dimensions (3D)], 2) topology (simply- vs multiply-connected), and 
3) shape (convex vs concave) of the FS. In particular, the term is absent for any quadratic (but 
not necessarily isotropic) spectrum both in 2D and 3D. The term is also absent for a convex 
and simply-connected but otherwise arbitrarily anisotropic FS in 2D. The origin of this nullification 
is approximate integrability of the electron motion on a 2D FS, where the energy and momentum 
conservation laws do not allow for current relaxation to leading -second-order in T/Ef {Ef is the 
Fermi energ y). If the ter m is nullified by the conservation law, the first non-zero term behaves 
as r^. The same applies to a quantum-critical metal in the vicinity of a Pomeranchuk instability, 
with a proviso that the leading (first non-zero) term in the resistivity scales as T 3 (T 3 ). We 
discuss a number of situations when integrability is weakly broken, e.g., by inter-plane hopping in 
a quasi-2D metal or by warping of the FS as in the surface states of topological insulators of the 
Bi2Te3 family. The paper is intended to be self-contained and pedagogical; review of the existing 
results is included along with the original ones wherever deemed necessary for completeness. 



I. INTRODUCTION 

A scaling of the resistivity vi^ith temperature (T) is 
considered as an archetypal signature of the Fermi-liquid 
(FL) behavior in metals. This result owes its origin to the 
Pauli exclusion principle which dictates that, at low tem- 
peratures, only those quasiparticles that reside within a 
width of order T near the Fermi energy participate in 
binary collisions. This argument, however, applies only 
to the inverse of the quasiparticle relaxation time 1 /tcc 
but not to the resistivity, p, per se: the scaling of 
the former does not necessarily imply that of the latter. 
A very simple example is a Galilean-invariant FL, where 
the electron-electron (ee) interaction does not affect the 
resistivity, although 1/too, as measured, e.g., by thermal 
conductivity, does scale as T^. The reason is that, since 
velocities of electrons are proportional to their respec- 
tive momenta, conservation of momentum automatically 
implies conservation of the electric current. In order to 
achieve a steady-state current under the effect of an ex- 
ternal electric field, a momentum relaxation mechanism 
is needed. 

Of course, the FL of electrons in a metal is not 
Galilean-invariant. In the presence of lattice, the cur- 
rent may be relaxed by Umklapp collisions,!^ which con- 
serve the quasimomentum up to a reciprocal lattice vec- 
tor: k -(- p = k' + p' + b. Umklapp processes are allowed. 
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however, if the incoming electron momenta k and p as 
well as the momentum transfer q = k — k' = p' — p are 
all of order b. These requirements are satisfied 1) if the 
Fermi surface is large enough, e.g., at least quarter-filled 
in the tight-binding case,'2iand 2) if the interaction is suf- 
ficiently short-ranged. In conventional metals, these two 
conditions are easily met due to a large number of car- 
riers and effective screening of the Coulomb interaction; 
thus Umklapp collisions occur at a rate comparable to 
1 /toc , and p (xT^. 

However, there arc situations when these conditions 
are not met; e.g., the first condition is violated in systems 
with low carrier concentration, such as degenerate semi- 
conductors, semimetals, surface states of three dimen- 
sional topological insulators, etc., and the second condi- 
tion is violated when a metal is tuned to the vicinity of 
a Pomeranchuk- type quantum phase transitioiP (QPT), 
e.g, a ferromagnetic QPT. A Pomeranchuk-type QPT is 
a g = instability of the ground state, manifested by a 
divergence of long-wavelength fluctuations of the order 
parameter. The effective radius of the interaction me- 
diated by the exchange of such fluctuations diverges at 
the QPT. One of the consequences of this divergence is 
the FL breakdown, as manifested by a non-Fermi-liquid 
(NFL) scaling I/tco c>c T'' with 7 < 1, but another one is 
the concurrent suppression of Umklapp processes. 

If Umklapps are suppressed (and the temperature is 
too low for the electron-phonon interaction to be effec- 
tive), current can be relaxed only via electron-impurity 
(ei) collisions. Still, the normal, i.e., momentum- 
conserving, ee collisions can affect the resistivity, if cer- 
tain conditions are met. The main purpose of this pa- 
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per is to summarize and analyze these conditions. The 
combined effect of normal ee and ei interactions does 
not necessarily lead to the dependence (or its NFL 
analog) of the resistivity. Whether this happens de- 
pends on three factors: 1) dimensionality [two dimensions 
(2D) vs three dimensions (3D)], 2) topology (simply vs 
multiply connected), and 3) shape (convex vs concave) 
of the Fermi surface (FS). The term is absent not 
only for a Galilean-invariant but, more generally, for an 
isotropic FL with a non-parabolic spectrum, as well as 
for anisotropic but quadratic spectrum. In 2D, the condi- 
tions are more stringent. In addition to cases mentioned 
above, the term is absent for a simply-connected and 
convex but otherwise arbitrarily anisotropic FS. The rea- 
son behind this is that the term arises from electrons 
confined to move along the FS contour such that, for 
the convex case, momentum and energy conservations are 
similar to the ID case, where no relaxation is possible. 

The issue of an interplay between normal ee and ei in- 
teractions has a long history, and it is beyond the scope 
of this paper to give a comprehensive review of the exist- 
ing literature; some aspects relevant to 3D metals are re- 
viewed in Ref.m Very briefly, the first notion that normal 
processes can affect the resistivity even in a single-band 
metal probably goes back to the paper by Debye and 
Conwell.l^ There is also a large body of work on normal 
collisions in multi-band metals, following the original pa- 
per by Baber,^ both at the phenomenological (reviewed 
thoroughly in Ref. 7) and microscopic* levels. That mo- 
mentum relaxation occurs differently in 2D as compared 
to 3D was pointed by Gurzhi, Kopeliovich, and Rutke- 
vich, first for the electron-phonoiP and then for the ee^ 
interactions. Maebashi and FukuyamsP^I analyzed an in- 
terplay between normal and Umklapp collisions for an 
anisotropic 2D FS and found that the normal collisions 
do not give rise to a term as long as the FS is con- 
vex. Rosch and Howell (Ref. [Hi) and Rosch (Ref. HSb) 
showed that a similar nullification happens for the 
term in the optical conductivity in a disorder-free 2D 
system. Chubukov and two of us (D. L. M. and V. I. 
Y.) generalized the analysis for a NFL near the Pomer- 
anchuk QPT-f^^l Scaling of the resistivity near a convex- 
to-concave transition was studied in Ref. '14', This paper 
expands on our recent worka^^^ and provides some more 
details. 

It is worth noting that the effects studied in this paper 
occur already within the semiclassical theory of trans- 
port that neglects quantum interference between ee and 
ei scatterings. Whether semiclassical description makes 
sense is one of the issues analyzed in the paper (cf. 
Sec. VA|: as a general rule, semiclassical effects can be 



considered separately from quantum-interference ones in 
the ballistic but not in the diffusive limit. Another ef- 
fect not captured by the semiclassical Boltzmann equa- 
tion (BE) is the viscous correction to the resistivity, dis- 
cussed in Sec. |VB[ The viscous correction is expected 
to be the leading contribution to the resistivity resulting 
from the ee interaction in the high-temperature, hydrody- 



namic regime, when the ee mean free path is sufficiently 
short. 

The rest of the paper is organized as follows. We begin 
by formulating the problem in terms of the Boltzmann 
equation (BE) both for the FL and the NFL cases in 
Sec. [n] In Sec. |III[ we solve the BE perturbatively with 
respect to ee scattering, which is an adequate approxi- 
mation at low enough temperatures, and analyze various 
stituations mentioned above. In Sec. IV we discuss the 
opposite limit of high temperatures, when the ee contri- 
bution to the resistivity saturates, and show that a true 
scaling regime, with an appreciable difference between 
the low and high temperature limits of the resistivity, 
does not exist in a single-band metal (Sec. IV A). Such a 
regime is shown to exist for a two-band metal with very 
different masses (Sec. IV B). In Sec. |vj we analyze the 
limits of the validity of the results based on the semi- 



classical BE with respect to both quantum (Sec. VA) 
and classical (Sec. VB| correlations between ee and ei 
interactions. 
SedVH 



Our concluding remarks are presented in 



II. BOLTZMANN EQUATION: GENERALITIES 
A. Collision integral 

The most straightforward way to find the effect of the 
ee interaction on the conductivity in the semi-classical 
regime is via the Boltzmann equation (BE) which, for 
the case of a time-independent and spatially-uniform ex- 
ternal electric field E, reads 



eE 



-loi [/k] — ^cc [/k] 



(2.1) 



where — e is the electron charge and /k is the distribu- 
tion function. The collision integrals Ice and Id on the 
right-hand side describe the effects of the ee and ei in- 
teractions, respectively. Explicitly, 



^oi = J Wk'k (/k - (Ck - £k') : 



(2.2) 



and 



W^k,p^k'p' 

IpJp' Jk' 

xS (ek + Ep - Ek' - Ep') ^ (k + P - k' - p') 

X Lfk/p (1 - ./k') (1 - /p') - /k'/p' (1 - h) (1 - /p)] , 

(2.3) 

where is a short-hand notation for J j^jdj £^nd iOk,k' 
and VFk.p-i-k'p' are the ei and ee scattering probabilities, 
correspondingly. For a weak electric field, the left-hand 
side of the BE reduces to evk • EnJ^., where Vk is the 
electron group velocity and nk = 7T,(ek) is the equilibrium 
distribution function, with prime denoting a derivative 
with respect to the electron energy, £k (measured from 
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the Fermi energy) . Linearizing the ee collision integral on 
the right-hand side with respect to the non-equilibrium 
correction to rik, defined as 

/k = - rn^ffk = ?^k + n-k (1 - '^k) ffk, (2.4) 
one obtain^ 



Ice — / Wk,p-i.k'p' 

J p J p' 'f 

X (ffk + 9p - 9k' - Sp') "k"-p (1 - ?^k') (1 - «p') 
X(5 (k + p - k' - p') 5 (ck + Ep - ek' - £p') •(2-5) 



B. Pomeranchuk quantum criticality 

In addition to the case of a generic FL, we will be 
also interested in a special but widely studied case of 
a FL near a Pomeranchuk-type QFT^ which breaks 
the rotational symmetry and/or topology of the FS but 
leaves the translational symmetries intact. Examples 
of such a QPT include ferromagnetic and electronic ne- 
matic transitions.'^^ As opposed to, e.g., charge-density 
waves and antiferromagnets, both the ordered and dis- 
ordered phases are spatially uniform, and the transition 
is manifested via the divergence of certain susceptibility 
at g = 0. Therefore, critical fluctuations near the QPT 
are long-ranged, and the effective interaction among elec- 
trons, mediated by these fluctuations, is of a long range 
as well. Since a FL is, in general, unstable with respect to 
long-range interaction, the quantum-critical region of the 
phase diagram near the QPT is characterized by mani- 
festly non-Fermi liquid (NFL) properties, such as a di- 
vergence of the specific heat coefficient. However, the 
long-range nature of the effective interaction has another 
aspect; namely, small-angle scattering at critical fluctu- 
ations effectively prohibits Umklapp processes which, in 
the absence of disorder, are necessary to render the re- 
sistivity finite. Therefore, the ee contribution to the 
resistivity can result only from an interplay between ei 
and normal ee collisions. 

In this Section, we briefly summarize the properties of 
the simplest model describing a QPT of the Pomeranchuk 
type: the Hertz-Moriya-MiUis (HMM) model.IiSl In this 
model, electrons are assumed to interact via an effective 
potential proportional to the divergent susceptibility of 
the order parameter. Details of the effective interaction 
depend on whether instability occurs in the charge or 
spin channel but, for our purposes, it suffices to model 
the interaction by a scalar function 



(2.6) 



where lyp is the density of states, (5 > is the "distance" 
to the critical point along the axis of the control param- 
eter (pressure, doping, etc.), and a i s the radius of inter- 
action in the critical channel. Since (2.6) can be derived. 



strictly speaking, only in the random-phase approxima- 
tion, one needs to require that kpa^ I (Ref.ilTj). [Alter- 
natively, one can assume that the coupling between elec- 
trons and critical fluctuations is weak;^^ results of these 
two approaches differ only by re-deflnition of parameters.] 
The imaginary part of UcS results from Landau damping 
of critical fluctuations by itinerant electrons. The corre- 
lation length of critical fluctuations ^ = a/VS diverges 
at the QPT, where S ^ 0. 



For the interaction in Eq. (2.6), the inverse quasipar- 



tide lifetime (the imaginary part of the self-energy) be- 
haves as I/tcc oc for T ^ TpL = vpa'^/^^ and as 
1/Tee oc T^/^ fo^. y ^ y^,^ ^ rpj^^ energy scale Tfl 
separates the FL and NFL regions of the phase dia- 
gram. Since the momentum transfers are small in both 
regions (q ~ for T <C Tfl and q ^ {uj/vpa'^)^^^ ~ 
{T /vfo^Y^^ for T <C Ifl), the transport scattering time 

(1/Tec) (g/fcj.)^ In 



the NFL region, cx T(^+2)/3- in the FL regime. 



is longer than the lifetime: l/r^, 

1 /tII cx with a small prefactor.SSl The conventional 
wisdom was that "transportization" of the relaxation 
time was the only manifestation of the long-range nature 
of the interaction, so that the resistivity could simply 
be obtained by substituting the transport time into the 
Drude formula. r^j^-g y^^^As p cx T^/^ and p cx T^/^ in 
the NFL regions in 3D and 2D, correspondingly. The 5/3 
scaling of the resistivity is indeed close to what has been 
observed experimentally in a number of itinerant ferro- 
magnets near a QPT.^il The reasoning tacitly assumes, 
however, that Umklapp collisions are still present, and 
occur at a rate comparable to that for normal ones, so 
that the transport time for normal collisions gives a rea- 
sonable estimate for the Umklapp scattering time. As we 
have already pointed out at the beginning of this Section, 
this assumption is not satisfied for a long-range interac- 
tion. In the next Section, we will quantify this statement. 
Before we proceed further, some general comments on the 
HMM model are in order. 

First, we are going to use the BE even in the NFL 
region of the QPT, where quasiparticles are not well de- 
fined, i.e., when I/tqc ^ T. This seems to be incon- 
sistent with the general criterion of the validity of the 
BE.I^ However, well-defined quasiparticles are not re- 
quired for the BE to be valid in a special case, when 
the effective ee interaction can be treated in the Migdal- 
Eliashberg approximation, i.e., when the self-energy de- 
pends on the electron energy but not the momentum and 
vertex corrections are small. In this case, the BE can 
be derived in the Keldysh technique without any condi- 
tions on the parameter t^cT , as long as a much weaker 
condition Teoe_F >• 1 is satisfied.t^^l This argument, for- 
mulated first by Prange and Kadanoff for the electron- 
phonon interaction,^'* was used later by a number of re- 
searchers in a wider context ,t251 and also applicable to 
NFL systems, given that they allow Migdal-Eliashberg 
description. Having said that, we come to the second 
point, which is that the Migdal-Eliashberg treatment of 
the HMM model, thought previously to be controllable 
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a) b) 




FIG. 1. (a) Umklapp process for a long-range electron- 
electron interaction. One of the electrons (with initial mo- 
mentum k) is scattered by a small angle via the ee interaction, 
while another one (with initial momentum p) is scattered by 
the lattice all the way across the Brillouin zone. 6) Umklapp 
processes for large momentum transfer. The original FS is in 
the center. Reprinted from Ref. 1131 courtesy of the APS. 



at least in the 1 /N approximationfiSl has recently been 
shown to break down beyond the second loop order.l^ 
While acknowledging this problem, we remark that the 
processes responsible for this breakdown are effectively 



ID-Uke scattering events, in which both the initial and 
final fermions move along the same line. Although these 
processes are dangerous for the single-particle self-energy, 
their contribution to the conductivity should be reduced 
by at least the "transport factor" , which discriminates 
against small-angle scattering. 

Once the BE is adopted, the difference between the FL 
and NFL regimes becomes formal: the dependence of the 
ee scattering probability on the energy transfer may be 
neglected in the former but not in the latter. 

C. Matrix elements on a lattice: normal vs 
Umklapp processes 

The interaction potential between electrons on a lat- 
tice [/(ri,r2) depends on the coordinates of two elec- 
trons separately rather than on their relative coordinate; 
transforming to the center-of-mass and relative coordi- 
nates, U (ri,r2) — U [ri — Y2, '"^'^'"^ ), where [/ is a pe- 
riodic function of ^^^^^ but not of Vi ~ T2- (The time- 
dependence of the effective interaction is not essential for 
the analysis below, and will be omitted.) Consequently, 

?7(ri,r2)= j ^e-'i-(''i"""=^)e-*'"-(''i+""^)/2[/(q,b), 
•^q b 

(2.7) 

where b is the reciprocal lattice vector and the volume 
of the system is put to unity. The matrix element of U 
on the Bloch wave functions (r) ~ X^b (b) e'^''"*')''" 
reads 



Mkp^k'p' = X! '5k-k'+bi-b2,p'-p+b3-b4"k' (bi) Uk (b2) Mp' (bg) Up (b4) C/ (k - k' + bi - b2 - b/2, b) .(2.8) 

b,bi ...b4 

Now we consider a long-range interaction, relevant for a FL near the Pomeranchuk instability. In this case, the matrix 
element is non-negligible only if the first argument of U is as small as possible, which means that b = 2(b2 — bi) and 
k w k': 

A4p^k'p' = ^k-k'+bi-b.,p'-p+b3-b4<' (bi) Uk (b2) u*^, (bg) Up (b4) [/ (k - k', 2(b2 - bi)) . (2.9) 

bi...b4 



We see that the condition for an Umklapp process be- 
comes very stringent.'^ since k « k', the momentum 
conservation condition k — k' = p' — p + b can only 
be satisfied at special points, where p' — p ~ b and 
b = b4 — bg -I- bi — b2 is just another reciprocal lattice 
vector. As Fig. ^p, shows, this is only possible if p and 
p' are located at the edges of the Brillouin zone (and 
the FS is open). The volumes (areas) around the spe- 
cial points are small-in proportion to a small momentum 
transfer q. The corresponding scattering rate is smaller 
than the transport rate of ee collisions by a factor of . 
For HMM criticality, where q oc T^/^, this implies that 
the contribution to the resistivity from the process de- 



picted in Fig. [T}j scales as T^t^+^'/s^ j g., as T^/^ in 
3D and as in 2D. In both cases, the exponents are 
larger or equal than 2. This means that the NFL con- 
tribution to the resistivity is smaller (3D) or comparable 
(2D) to the FL (T^) contribution, arising from Umklapp 
scattering in the channels that are not affected by the 
proximity to a QCP, e.g., from the charge channel in the 
vicinity of a magnetic instability. In addition, processes 
in Fig. are, in fact, " pseudo-Umklapps" because they 
can be viewed as normal processes on a closed (hole) FS. 
The "real" Umklapps, shown in Fig. [T^, can occur only 
if the constraint of small momentum transfer is relaxed 
or else, near half-filling, when the "gap" between the FS 
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and the edges of the Brillouin zone is small. Half-filling, 
however, is more likely to result in a finite-g instability 
of the ground state, e.g., antiferromagnetism, rather than 
in a Pomeranchuk QTP. From now on, our analysis will 
be focused on normal ee collisions, the matrix clement of 



which is given by Eq. (2.9) with bi — b2 = — b4. The 



corresponding scattering probability, averaged over spins 
of the initial states and summed over spins of the final 
states, reads 



W, 



kp— >k'p' 



•k'p' 



►p'k'l 



kp— J^k'p' I 



kp— »p'k' I 



-Re(Mkp^k'p'M*p^p,k- 



(2.10) 

Wkp-fk'p' has certain symmetries. First, we assume the 
microreversibility property 

W^kp^k'p' - VKk'p'^kp. (2.11) 
In addition, since electrons are indistinguishable, 

W^k,p-i-k'p' = W^kp^p'k' = W^pk^k'.p' = W^pk^p'k'- 

(2.12) 



Finally, combining (2.111 and (2.121, we obtain 



W^kp-*k'p' — VKk'p'^kp — Wp'k'^pk 



(2.13) 



With only normal ee collisions taken into account, the 
total electron momentum is conserved, i.e.. 



0. 



(2.14) 



Notice that although /ei is written down in its most 
gene ral form that holds true as long as Wk,k' obeys uni- 
taritjJ^^'^IISSl^ in writing down /ee we have already as- 
sumed that W^k,p->k'p' obeys the microreversibility con- 
dition \2.ll\. 



D. General properties of the solution 

Before proceeding with a more detailed analysis of the 
BE, we make a few general comments. 

1. "Hidden" phonons. The linearized form of 
the steady-state BE assumes implicitly that the 
electron-phonon interaction is also present in the 
system; otherwise, the total electron energy will 
increase indefinitely due to the work done by the 
electric field. As usual (see, e.g., Ref. [7|, we assume 
that the temperature is low enough so that one can 
neglect a direct electron-phonon contribution to 
the resistivity (which requires that Tcph S> t^c , Tgi , 
where r's are the transport scattering times for cor- 
responding processes) but high enough so that, for 
a fixed electric field, the electron-phonon interac- 
tion can still equalize the electron and lattice tem- 
peratures (which requires that the work done by 
the electric field on the energy relaxation length is 
much smaller than the temperature). 



2. Parity of a non-equilibrium part of the dis- 
tribution function. A linear in E term in /k can 
be written as 



<^/k = /k - f^k 



k' 



(2.15) 



where Ak contains explicitly only the effects of the 
ee and ei interactions. At low enough tempera- 
tures (as specified in the previous paragraph), the 
electron-phonon interaction shows up only in the 
next-quadratic-term and does not affect the resis- 
tivity directly. At even lower temperatures, when 
Tec ^ Ti, ee scattering can be treated as a perturba- 
tion to ei scattering. In this case, Ak is determined 
by the crystal symmetry and by the ei scattering 
probability and, without specifying both of them, 
no properties of Ak can be further inferred. How- 
ever, if the ei scattering probability satisfies the mi- 
croreversibility condition,^ i.e., Wk,k' — Wk'.k, then 
Ak is odd in k. Indeed, reversing the sign of k in 
the BE and relabeling k' — > — k', we obtain 



evk • E^k = 



-k'.-k (/-k' - ./-k) 5 (ek - Ek') ■ 



(2.16) 

Using time-reversal symmetry i«k,k' — ^^-k',-k 
(which is guaranteed in the absence of the mag- 
netic field and magnetic order) and microreversibil- 
ity, we see that A_k — — Ak. This is the prop- 
erty of the non-equilibrium distribution function 
we will be using later on. To simplify the presenta- 
tion, we will first use a model form of the ei colli- 
sion integral, namely, a relaxation-time approxima- 
tion /ei = (/k — ^k) /ti, (cf. comment in Ref. 'M'), 
which allows for a closed-form solution, and then 
extend the proof for the general form of /k given 
by Eq. ( pls] ). 

However, one has to keep in mind that, beyond 
the Born approximation, microreversibility is not a 
general principle but a consequence of two micro- 
scopic symmetries, i.e., symmetries with respect to 
time- and space-inversions, and is thus absent in 
non-centrosymmetric systems.l^ 

No disorder no steady-state linear-response 
regime. Since the momentum is conserved in nor- 



mal collisions, the collision integral (2.5 1 is nullified 
by a combination B • k, where B is k-independent 
but otherwise arbitrary. This means that there 
is no unique steady-state solution in the linear- 
response regime. Obviously, the steady-state solu- 
tion is absent because the total momentum of the 
electron system (per unit volume), K = Jj^k/k, in- 
creases with time. Indeed, restoring the time and 
spatial derivatives in the BE, multiplying it by k 
and integrating over k we obtain 



dt 



9H,. 
dx4 



(2.17) 
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where Hij = Jj^kiVjfk. Integrating by parts in the 
right-hand side and taking into account that the 
number density N = J^. /k, we obtain 



dt 



an,. 

dxi 



-eNE, 



(2.18) 



The left-hand side is just the continuity equation 
while the right-hand side is the total force per unit 
volume. Therefore, although the electron liquid is 
not, generally speaking, Galilean-invariant, it is ac- 
celerated as a whole by the electric field. (In a crys- 
tal, an increase of the momentum in time leads to 
Bloch oscillations of the current; the current aver- 
aged over time is equal to zero.) Therefore, one 
needs to invoke impurity scattering in order to ren- 
der the problem well-defined.'^ 

No lattice no term in the resistivity. 

Adding just disorder but no lattice does not give 
rise to a term in the resistivity. Notice that this 
statement is weaker than "the ee interaction does 
not effect the resistivity at all", which is true if 
Wk'.k depends only on the scattering angle but not 
on the electron energy. The simplest case is that 
of point-like impurities, when Wk.k' = l/i^(ek)Ti, 
where v{si^) is the density of states (per one spin 
component) and ri is a constant. In this case, the 
BE reduces to 



evk • EnL = - 



./k - ./ 



(2.19) 



where 9 {9') is the angle between E and k (k'), and 
Wo is an odd function of the polar angles that van- 
ishes when substituted into the collision integral. 
In terms of angular harmonics, the BE reduces to 



ewk 



-V{cos 19) - /oc(2.23) 



with 



(£k). 



(2.24) 



If Tj^^^(£k) does not depend on the electron en- 
ergy, one proceeds in the same way as for point- 
like impurities, i.e., one obtains a direct relation 
between j and E by multiplying (|2.23 1 by v and 



integrating over k. The resulting (Drude) conduc- 
tivity a — ne^rf^ /m contains the transport time 
r*"" = T^^^ but no effects of the ee interaction 
(again, up to FL renormalizations) . 

If w does depend on the electron energy, as it is of- 
ten the case for semiconductors, the integral k/oi 
does not reduce to the electric current, and one 
needs to solve for /k in order to find the conduc- 
tivity. Since the ee interaction affects /k, the con- 
duct ivity is also affected. However, as we show in 



Sec. Ill C 3 the effect leads only to a T term in 



the resistivity (or InT in 2D). 



where / is an average of /k the directions of k. 
In the absence of lattice, k = mv and hence the 
electric current j = — 2e /j^v/k — — (2e/m) /j^k/k. 
Now one can multiply the BE equation by v and 
integrate over k, upon which Joo-in accord with 
(2.14)- drops out, and obtain a relation between j 



and E directly, without solving for /k: 



f . m 

e / Vk Vk • E rik = - — J 
Jk 2eTi 



(2.20) 



The resuling conductivity a = ne^Ti/m does not 
contain any effects of the ee interaction, except for 
FL renormalizations of m and Tj. 

The same is true if the scattering probability de- 
pends only on the angle between k and k'. Param- 
eterizing Wk.k' as 

Wk^k' =w(£k,k-k')/i^(ek), (2.21) 

we expand /k and iu ^£kik- k'^ over a complete 
basis set of, e.g., Legendre polynomials in 3D: 



/k 



= ^ w^^^ (ek)'P^ (cos 9)Ve{cos 9') + Wo, (2.22) 



III. ELECTRON-ELECTRON CONTRIBUTION 
TO THE RESISTIVITY 

A. Do normal ee collisions affect the resistivity? 

It may seem that the reverse statement to the head- 
ing of item ^4 in the previous Section ("no lattice-no 
term in the resistivity" ) should be "a term in the 
resistivity occurs in the presence of both disorder and 
lattice". Indeed, while disorder takes care of momentum 
relaxation, lattice breaks the Galilean invariance. As a 
result, Vk = dei^/d\i ^ k/m, which means that momen- 
tum conservation does not imply current conservation, 
and one cannot obtain a relation between the current 
and the electric field without actually solving the BE. In 
general, therefore, one should expect a term in the 
resistivity. While it is really the case in 3D, it turns out 
that the conservation laws in 2D forbid the term for 
a convex and simply-connected but otherwise arbitrary 
FS. 



B. Low temperatures: perturbation theory 

In this Section, we discuss the case of low tempera- 
tures, when the ee collisions are less frequent than the ei 
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ones. In this case, the ee contribution can be found via 
the perturbation theory with respect to Ice- We begin 
with the simplest-isotropic- model for electron-impurity 
scattering, when the BE is given by (2.19). However, 



we keep the dependence of the ei relaxation time on the 
electron energy for the time being. 



yields 



5k^ =eTi(£k)vk-E/T. 



(3.1) 



Next, we substitute g^^ back into (2.19) and find a cor- 
rection due to Icr 



(2) _ n (£k) J 

" Tn', ' 



9^' 



(3.2) 



At the first step, we solve (2.19) with Joe = 0, which 



The corresponding correction to the ij component of the 
conductivity tensor is given by 



5 a. 



T 



dujdeit^der 



^^^M^k.p(q,c.)^LA^^ 



xn (£k) n (ep) [1 - n {e^ - w)] [1 - n {Sp + uj)] 6 {e^ ~ £k-q - w) (5 {Sp - e. 



p+q 



■ w) 



(3.3) 



Here, q = k — k' = p' — p is the momentum trans- 
fer, dflk is the surface (line) element of an isoenergetic 
surface (contour) at energy in 3D (2D), and A£ = 
Ti(£k)vk + Ti(ep)vp - ri(ek - w)vk-q - Ti(£p + a;)vp+q 
is a vector measuring the change in the total "vec- 
tor mean free path " ^k = VkTi(ek) due to ee colli- 
sions. The energy transfer was introduced by re- writing 



the energy conservation law as 5 (ck + Sp — Sk' — £p' ) ~ 
J dioS (ck — £k' — i^) S {Sp ~ Ep' + Lu). The scattering 
probability P1/k,p(q, i^) = W^k,p->-k-q,p+q is now allowed 
to depend on uj. Using the symmetry properties of 
Wk.p(q, w), one can cast (3.3) into a more symmetric 
form 



San = 7^ 



d^q 



diodekdsr 



dflk da 

Vk Vp 



''Wk,p{q,u;)AeA£^ 



xn (ek) "-P (cp) [l-n {e^ - w)] [1 - n {sp + lo)] S (sk - £k-q - uj) S {Sp - £p+q + uj) . 



(3.4) 



Now let us count the powers of T in (3.4). Each of the 



three energy integrals (over uj, e^, and Sp) gives a factor 
of T which, in a combination with the overall 1/T factor, 
already gives a dependence, as is to be expected for 
a PL. The result holds as long as the integral over q 
does not introduce additional T dependence. This is the 
case in the FL regime, when typical q are of order of the 
ultraviolet cutoff of the problem, i.e., the smallest of the 
three quantities: the reciprocal lattice vector, a typical 
size of the FS, and the inverse radius of the ee interac- 
tion. In this case, the uj dependence of Wk,p (q, uj) can 
be neglected. The energy dependence of t; contributes 
only to higher order terms in T and we neglect it for the 
time being as well, so that A£ — TiAv with 



Av = Vk + Vp - Vk-c 



(3.5) 



being the change in the electron current due to ee colli- 
sions. Finally, since the integrals of the combinaton of 
the Fermi functions over energies already produce a fac- 
tor of T^, electrons can be projected onto the FS in the 
rest of the formula. This means that one can drop uj in 



both (5-functions and perform the surface integrals over 
the FS. 

We pause here to remark that neglecting w in the 
(5-functions does not mean performing an expansion in 
w/ck, uj/sp, etc. In fact, all quasiparticles energies (ck, 
Ek-q, etc.) are equal to zero because the electrons were 
projected onto the FS! What it really means is that the 
(5-functions impose constraints on the angles between k 
and q (and p and q) with electrons' momenta being on 
the FS. Typical values of these angles are determined by 
the ratio of typical q {= q) to kp. In a system with a 
short-range interaction, q ~ min{fc^, 1/ao}, where oq is 
the lattice spacing; therefore, typical angles are of order 
unity. On the other hand, typical uj (= uj) are of order T, 
and corrections to angles due to finite uj are small as long 
as T ^ mm{ep, W}, where W is the bandwidth; the last 
condition is implied anyhow to be in the FL regime. If 
the interaction radius, tq, is much longer than both the 
lattice spacing and the Fermi wavelength, q is small but 
in proportion to l/rp rather than to T, while <D is still 
of order T. This means that effective ultraviolet energy 
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scale is reduced to vp/ro, and the FL description is valid 
only at low energies, where the effect of a finite energy 
transfer on the kinematics of collisions is negligible. This 
can be illustrated for a simple example of the quadratic 
spectrum, when the angle between, e.g., k and q, satis- 
fies cos 0k, q — [q^ /2m + uj)/vpq. Neglecting w is justified 
as long as r ^ (f /2m. Notice that this simplification 



ductivity in all dimensions _D > 1, because the small q 



is valid even near QPT (cf. Sec. II B), where the scaling 
dimensions of w and q are different: u) ^ T but q cx T^^^. 
A characteristic temperature, below which the condition 

satisfied, coincides with the 
scale below which the quasiparticle description breaks 
down, which is the regime of main interest for quantum 
phase transitions. 

As another remark, typical q may be different for dif- 
ferent observables. What we said above is true for the 
leading term in the ee contribution to the electrical con- 



behavior of the integrand in (3.4 1 is regularized by the 
A£^ factors that vanish in the limit of g — > 0. (This 
is analogous to the regularizing effect of the 1 — cos 9 
factor in a transport cross-section for elastic scattering.) 
However, when calculating the single-particle lifetime 
(the imaginary part of the self-energy/S2l and thermal 
conductivity'^^ in 2D, one runs into infrared logarithmic 
divergences, which means that the infrared region of the 
momentum transfers {q T/vp) does contribute to the 
result. In those cases, neglecting uj in the J-functions is 
not justified. (The subleading terms in the conductivity 
also require more care; see Sec. IIIC3 below.) 

Coming back to the main theme, we focus now on the 
FL case, when one can also neglect uj in the scattering 
probability. After all these simplifications, the diagonal 
component of the conductivity reduces to 



Saji = T: 

1 1 2^ 1 



d^q 
i2nf 



1 p J F 

"""^ l^k.p (q, 0) (Av^^,^ ^ (g^) _ ^ _ _ ^ 



C(5 (ek - Ek-q) |ek=0'^ ('^p - ^p+q) |ep=0, 



(3.6) 



r 



where superscript F indicates that the corresponding explicitly and, instead of Eq. (3.8), we obtain 
quantity is evaluated at the FS. Now the integrals over 
all energies can be performed with the help of an identity 



dex I de2 I dujn[ei)n{e2)[^ — n[ei — u) 



X [1 - n (£2 + ^)] = 



2ti' 



(3.7) 



d^q 
J2nf 



daedal 

P i?k,p (q) 



X [Avf J (ek - ek-q) l.^^o-^ (fp - Cp+q) |.p=o, (3.9) 
where 



and we obtain a term in the conductivity with a pref- 
actor given by a certain average over the FS 



d^q 
J2nf 



dal da^ 



|^l^k,p(q,0) 



X (Av*)2,5 (ck - Ek-q) \e^=o5 (£p - £p+q) |ep=0 



(3.8) 



i?k,p (q) = j dio {uj^T^) Wk,p (q, uj) N{uj) [N{uj) + 1] 

(3.10) 

and N{uj) is the Bose function. For the effective interac- 



tio n fro m Eq. (2.6 1 at the QPT (d = 0), power counting 
of (3.9) gives San cx T^^+^'Z'^, which coincides with the 



estimate based on the transport time (cf. Sec. II B). As 



in the FL case, however, one needs to make sure that the 
prefactor is non-zero. 



Clearly, whether the leading correction to the residual 
conductivity indeed scales as depends on whether the 
integral over the FS is nonzero. Since the integrand is 
positive, the integral may vanish only if Av = under 
the energy conservation constraints imposed by the S- 



functions. As a simple check, we apply Eq. (3.8) for the 
Galilean-invariant case, when Vk = k/m. In this case, 
Av = 0, as it should be. 

We now consider a more general situation, when the 
w dependence of the scattering probability is important, 
which is the case, e.g., near a QPT. In this case only 
two out of the three energy integrals can be performed 



C. Cases when the leading term vanishes 

1. Isotropic system with an arbitrary spectrum 

The first case is that of an isotropic but otherwise ar- 
bitrary energy spectrum. Such a situation may arise due 
to relativistic effects. Another (pseudo-relativistic) ex- 
ample is weakly doped graphene with a negligibly small 
trigonal warping of the FS. Since £k is a fimction of |k| 
only, the (5-function constraints in Eq. (3.8) imply that 
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k' = k 4 



k ^ — ' p 

Coop er channel velocity swapping no scattering 




p =p 



Av = 



Av = 



Av = 



FIG. 2. Isotropic case in 2D: three possible scattering pro- 
cesses none of which leads to current relaxation. 



|k - q| and |p| = |p + q|. Then, 



dsk 



■3 1 



iP — 2 
•'^'1 ~ 9(/fc2) 



X (k] -qj) ^^ik){kj ~qj), 



|k-q| = |k| 



(3.11) 



where ^(fc) = v-^/k. Notice that the second line in 
Eq. ( 3.11 ) is not an expansion in small q but an exact rela- 



tion. Substituting Eq. (3.11) (and similar expressions for 



and Wp+q) into Av, it is easy to see that Av vanishes 
identically. Thus, there is no correction to the resis- 
tivity of a non-Galilean-invariant but isotropic system. 
This result also holds for a general quadratic spectrum 
£k = kikj /2mij, in which case Vj = ki/rriji and Av-' = 0. 

Notice that, in contrast to the Galilean-invariant case 
(with Ek = k"^ /2m — ep), when not only the term but 
all higher order terms are absent, higher order (T"*, etc.) 
terms are non-zero for a non-parabolic spectrum. 



2. Approximate integrability: convex and simply connected 
Fermi surface m 2D 

Kinematics of ee collisions on a circular FS. 

The fact that term in the resistivity is absent for an 
isotropic FS does not mean that it is necessarily present 
for an anisotropic FS. In fact, the term is also absent 
for a simply-co nnected and convex but otherwise arbi- 
trary FS in 2D! ^° I ^^ I ^ . Before considering the general 
case, however, let us study the simplest example of such 
a FS, i.e., a 2D circular FS with quadratic spectrum. 
Since this is just a Galilean-invariant case, we already 
know (cf. Sec. II D) that the ee interaction has no ef- 
fect on the resistivity. However, it is instructive to see 
in a geometrical way how the term vanishes-this will 
be useful for the subsequent analysis of the general case 
in 2D. Geometrically, one needs to find two initial mo- 
menta, k and p, belonging to the FS, such that the final 
momenta, k — q and p + q, also belong to the FS. As 
shown in Fig. 2] only three situations are possible!^ 1) 
Cooper channel, when the total initial and, therefore, the 
total final momentum is equal to zero; 2) swapping of ve- 
locities, when the initial momentum of one the electrons 
coincides with the final momentum of another electron 
and vice versa, i.e., p = k — q; 3) no scattering- this is 
the trivial case where the initial and final momenta of 
individual electrons are the same. For all of these cases, 
Av = and thus the term is absent. To see that 
these situations indeed exhaust all the possibilities, one 
can solve the momentum and energy conservation equa- 
tions, i.e., k— k' = p' — p = q and k^ — k'^;p^ — p'^, sub- 
ject to the additional constraint k — p — k' — p' — kp- 
This leads to two equations: — 2kqcos9kq — and 
q^ + 2pqcos9pq = 0, where 9ij denotes the angle between 
the vectors i and j. The three possible solutions are: 
9kq — 9pq = TT, Corresponding to case 1); 0kq + dpq — tt, 
corresponding to case 2) ; and q = for arbitrary Okq and 
Opq , corresponding to case 3). 

Kinematics of ee collisions on a generic convex 
FS. The situation described above is not specific to a 
circular FS in 2D but occurs also for a generic convex FS, 
see Fig. [3^. Indeed, introducing a new variable p = — p 



and symmetries of the scattering probability, we obtain 



3.8 ) and using the time-reveral symmetry (e_ 



5^^^ = -y e^rV^ I ^(jicji^^w^,^ (q, 0) [vj, - vl - + vi_j' 5 {e^ e^-q) 5 (sp - ep_q) . 

(3.12) 

I 



For given q, we must find two momenta satisfying the 
relations Ek — £k-q and ep — £p-q- Geometrically, find- 
ing the solution to these two equations is equivalent to 
shifting the FS by q, and finding the points of intersec- 
tion between the original and the shifted FSs. A convex 
FS has at most two self-intersection points. Therefore, 



the equation £k = £k-q has only two solutions. In ad- 
dition, if k is a solution, then — k -I- q is also a solu- 
tion so that the roots of the first equation form a set 
{k, — k -I- q}. Since the second equation is the same, its 
two roots {p,— p + q} = {— P,P + q} must coincide with 
the roots of the first equation. This can happen if 1) 
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k = — p, which gives the Cooper channel or if 2) k = p+q 
which gives swapping. The situation with g = 0, when no 
scattering occurs, is trivially possible. For all the scat- 
tering processes listed above, Av — and the term 
vanishes. 

Being purely geometrical, the preceding analysis is 
equally valid for the NFL case, with the conclusion that 
the T^/^ term vanishes as well. 

Beyond the relaxation-time approximation. 
Although the analysis above was based on Eq. (3.121, 



obtained in the relaxation-time approximation for ei 
scattering, it can be readily extended for the general 
form of the ei collision integral in Eq. (2.2). The non- 



equilibirum part of the distribution function in the pres- 
ence of ei scatter ing only is given by (2.15), which im- 
plies that g^^ in (2.4) is replaced by g^^ = Ak ■ E. The 
lowest-order iteration in ee scattering is to be found from 
an integral equation 



1 



Tni 



/cc[Ak-E<]. 



(3.13) 



The ei collision integral can be viewed as an integral op- 
erator, the inverse of which is defined by 



k 



k,k'/k 



(3.14) 



system behaves similar to a ID system, where binary col- 
lisions do not lead to relaxation. The analogy works be- 
cause, to find the leading (T^) term in the conductivity, 
it suffices to project electrons onto the FS, which is a line 
in 2D. Therefore, kinematics effectively becomes ID and, 
although this is a 2D case, we have an integrable system. 
However, this analogy has certain limitations. First, the 
2D case is integrable only with respect to charge but not 
thermal current relaxation, whereas there is no relaxation 
of all physical quantities in ID. Second, even the charge 
current relaxation is absent only up to next-order-terms 
in T/ep (see Sec. |IIIC3| . Third, not any FS line in 2D 
is integrable: concave and multiply-connected contours 





FIG. 3. (a) A convex contour has at most two self-intersection 
points (marked by dots), (b) A concave contour can have 
more than two self-intersection points (six in the example 
shown) . 



where k = k/fc. Thanks to microreversibility, Ok.k' 
Ok'.k- A formal solution of (3.13) is 



(2) 

5k 



1 



Tn' 



-/Ti/,,[Ak-E]. 



(3.15) 



behave in a non-integrable way. With all these limita- 
tions in mind, we will refer to the 2D convex case as to 
" approximate integrability" . 



Using the microreversibility property of Ok,k' and the 

(2) 

fact that Ak is odd in k, it is easy to see that g^^ is odd 
in k as well. This is all one really needs to repeat the 
steps of the previous analysis. The correction to the con- 
ductivity now contains a combination Av'AA-', where 



AA 



Lk -f- Ap — Ak' — Ap/ . Being odd in all momenta, 
AA behaves in the same way as Av upon the change 
p — > —p. The scattering processes are classified in the 
same way as before, and the vanishing of the term 
follows from the vanishing of Av. 

Approximate integrability. A limited number of 
possible outcomes of the ee collisions means that our 2D 



Subleading corrections to the resistivity when the leading 
term is absent. 



Higher order term from electrons away from the 

FS. To find the subleading correction for the case con- 
sidered in the previous section, we go back to Eq. (3.4 1, 



replace again T{ by a constant in the scattering probabil- 
ity, but now, instead of neglecting uj in the 5 functions, 
expand the product of the 5 functions to second order in 
UJ. The zeroth-order term, 5 (sk — £k-q) <5 (£p — £p+q), 
nullifies Av'. The odd in lo terms vanish upon integra- 
tion over ek,£p, and uj. In the FL case, this gives 



5ai 



1^- 
2 T 



da\^ da-p 



M^k,p(q,0) 



X [Av*] ^ n (£k) n (ep) [1 - n (sk - t^)] [I ~ n {e^ + lo)] 

(ek - ^k-q) (£p - ep+q) - ^ {(5" (ek - ^k-q) S {sp - Cp+q) + (5 (ek - ^k-q) S" (ep 



'P+q 



)} 



(3.16) 



The derivatives of the i5-functions produce the same roots tegrating by parts, we make the derivatives to act on 
for k and p as the (5-functions themselves. However, in- 
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[Av*] ^ . Although [Av*] ^ vanishes for k and p satisfy- 
ing energy and momentum conservations, its derivatives 
do not. This makes the integral non-zero. Since we now 
have two more factors of cj the correction to the conduc- 
tivity scales as 



San oc T 



(3.17) 



In more detail, let ko be one of the roots of the equation 
Ek = Ek-q- The corresponding root for p is then po = 
kg — q. Expanding Av around the roots gives 



Av* = {[6k-Sp] ■ V) (vj 



ko 



'ko-q 



where (5k = k — kg and (5p = p — ko 
integration proceeds as in the integral 

1 



) , (3.18) 
Subsequent 



dxdyS' (x) 5' (y) ^ {x 



y) ——J dx6'{x)x — l, 
^ (3.19) 

where 6k. and dp play the roles of x and y (and similarly 
for an integral with a product 6" {. . .)6{. . .)). Further 
cancelations for a particular FS may make the T depen- 
dence even weaker but the generic answer is T^. 

Clearly, going away from the FS produces an extra fac- 
tor of . Since w ~ T in the NFL as well, the result for 
the NFL regime is obtained by multiplying the "naive" 
estimate San oc T^/a by T^, which gives San oc T^°/^. 
This is obviously subleading to the term resulting 
from the FL interaction in non-critical channels. 

Energy-dependent electron-impurity relaxation 
time. In addition to the mechanism described above, 
there are other sources of higher than corrections to 
the conductivity; one of them is the energy dependence 
of Ti which we have neglected so far. This mechanism 
operates even in a Galilean-invariant system: although 
ee collisions conserve the momentum, they redistribute 
electrons in the energy space and thus affect the conduc- 
tivity, if Tj depends on the energy.^^ To e stimate the mag- 
nitude of this effect, we apply Eq. (3.4) to the Galilean- 
invariant case (v = k/m) and expand the impurity re- 
laxation times entering the "vector mean free path" as 



n (£i) = n{0) + t/ei, 

yields 



where t- = 9ri(ei)/i9ei| 



=0- 



This 



A£ = ^ (k - p - 2q) , 
m 



(3.20) 



Since (3.3) contains two factors of Ai, and each of them 
is proportional to u, we have an extra factor in the 
integrand. In 3D, this immediately gives a term 



Sa-^ (X iT!fT\ 



(3.21) 



In 2D, the situation is more delicate because the part of 
the integrand associated with the k — p term in (3.20) 
is logarithmically divergent. This a well-known "2D log 
singularity" that occurs, on a more general level, as the 
mass-shell singularity of the self-energy (see Ref. [32] and 
references therein) . This is also the same singularity that 





FIG. 4. (a) A 3D FS has an infinite number of self-intersection 
points (a line), (b) A multiply connected FS has more than 
two self-intersection points. 



one encounters when calculating the thermal conductiv- 
ity in 2D (in the absence of impurity scattering). '^'^ In- 
deed, our problem bears a formal similarity to that of the 
thermal conductivity because the change in the thermal 
Vk£k due to ee collisions 



current 



Jk 



Jp 



Jk 



Jp+q 



(k - p - 2q) + q(£k - gp) 

m 

(3.22) 



contains the same term as A£ in (3.20). The singularity 
can be resolved by the same method as in Ref.l33l i.e., by 
considering a dynamically screened Coulomb interaction. 
The result is that, similar to the thermal conductivity, 
the conductivity contains an extra log factor as compared 
to the 3D case: 



'2D 



o,{r[fTHu{ep/T). 



(3.23) 



The "2D log" does not occur in the term in the con- 
ductivity, if the latter is finite due to broken integrability, 
which is the subject of the next section. 

An extension to the NFL case is again, trivial, and we 
will not repeat the argument here. 



D. Non-integrable cases 

Concave FS in 2D. It follows from the previous 
discussion that whether the term is absent or present 
depends entirely on the FS having two or more than two 
self-intersection points. A concave FS in 2D can have 
more than two self-intersection points (cf. Fig. |3]d), 
therefore there are more than two solutions for the initial 
momenta for given q. Some of these solutions still corre- 
spond to "integrable" processes, encountered already for 
a convex FS, but the remaining ones do relax the current. 
Therefore, a term survives in this case. 

3D FS. In 3D, the manifold of intersection between 
the original and shifted FSs is a line, see Fig. |4][a). 
Therefore, the equation = £k-q has infinitely many 
roots. There is no correlation between the roots of the 
equations Sk = ek-q and £p = £p+q- Geometrically, this 
means that the initial momenta, k and p, do not have 
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to be in the same plane as the final ones, k' and p'. 
Therefore, an anisotropic (but not quadratic) FS in 3D 
allows for a correction to the resistivity. 

The T^/-^ term in the NFL regime survives for the same 
reason as well. Therefore, our theory at least does not 
contradict the experiments^ where such term was ob- 
served. 

Multiply connected FS. If the FS is multiply con- 
nected, a term in the resistivity is present, even if 
the individual FS sheets do not allow for a term on 
their own. Even more so, the individual sheets can even 
be isotropic. The reason is obvious from Fig. |4][b) which 
shows an example of two circular FSs in 2D. Clearly, the 
equation = £k-q has more than two roots even in this 
case. Thus, according to our previous arguments, there 
is no general reason for the vanishing of the term in 



into the in- and out-of-plane parts as 



such a situation. In Sec.|IV|and Appendix [A| we discuss 
the two-band case in 2D in more detail. 



E. Weakly-integrable cases 

In this section, we consider two situations when inte- 
jrability is broken only weakly. 



1. Quasi-2D metal 

The first case is a layered metal with a quasi-2D spec- 
trum which, for simplicity, we assume to be separable 



~-kz ' 



(3.24) 



where k|| and kz are the in-plane and out-of-plane com- 
ponents of the momentum, correspondingly. In the tight- 
binding model with nearest-neighbor hopping, = 
t±^ [1 — cos (fc^c)], where c is the lattice spacing in the 
z-direction. The metal is in a quasi-2D regime when 
tj^ <C Ef- In regard to the in-plane part of the spec- 
trum, e[| , we assume that the corresponding energy 
contours are anisotropic but convex so that, in the ab- 
sence of intcr-planc hopping, the T^-term in the in-plane 
conductivity would be absent. (If the planes are as- 
sumed to be Galilean- invariant, i.e., e^^^ ~ fcj^/2r7i||, as 

in a "corrugated cylinder model" , the T^-term is triv- 
ially zero because the in- and out-of-plane components of 
the momentum are conserved independently, and hence 



f^l - v", = 0.) To find the T^-term in the 

K,, P|| 



in-plane conductivity, we use a method similar to that 
in Sec. |III C 3[ i.e., we expand the (S-functions, except for 
that now we expand both in oj and ef, . As we explained 
in Sec. [Ill Bj the expansion in lo is really an expansion in 
Lj normalized by the appropriate ultraviolet energy scale 
of the problem. Likewise, the expansion in is really 
an expansion in t^/tp, which is a natural small parame- 
ter for a quasi-2D system. The zeroth-order term {lo = 0, 
£^ — 0) nullifies Av". The first-order terms also van- 
ish: the ones, proportional to oj, do so by parity, and the 
ones, proportional to £| , do so because the first-order 
derivatives of the 5-functions nullify (Av")^ after a sin- 
gle integration by parts. Finally, the cross products in 
second-order terms, being odd in oj, also vanish. There- 
fore, the only surviving second-order term is 
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(3.25) 



Equation (3.25) contains two independent corrections. 
All terms proportional to u'^ produce a correction 
to the conductivity that exists even in a purely 2D sys- 
tem. All terms containing the squares of the out-of-plane 
dispersions produce a correction.liSl Therefore, 



6cri, 



AaT^ 



(3.26) 



where i = x,y, and constants A4 and A2 depend on 
details of the in-plane spectrum; generically, A4 ~ A2. 
Equation (3.261 describes a dimensional crossover from 
the 2D-like regime {Sa^i cx T^) at T :$> t±_ to the 3D-like 



regime {San cx T^) for T < tj_). Notice that, in the 3D 
regime, the T^-term in the in-plane conductivity depends 
on the out-of-plane hopping. 



In the NFL regime, Eq. (3.261 is replaced by 
Sau = AiT^"/^ + A2tlT^/'. 



(3.27) 



2. Conductivity near the convex-concave transition 

In this Section, we consider a FS near a convex- 
concave transition which occurs when the Fermi en- 
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FIG. 5. (color online), (a) Even in the concave case, there 
can be only two intersection points if q is not along a special 
direction, (b) For q larger than a critical value, there are only 
two intersection points. 



ergy goes above a certain threshold value Such 
a situation is encountered, e.g., in the case of surface 
states of the Bi2Te3 family of 3D topological insulators, 
where the electron spectrum can be approximated by2SI 
= ±y^vj,k'^ + A2fc6cos2(36'), with 9 being the polar 
angle. Corresponding isoenergetic contours are shown in 
Fig. |6ja). We will be interested in the vicinity of the 
convex/concave transition, when je^ — £c| = |A| ^ Sc- 

We first consider the case of A ^ T, when the isoen- 
ergetic contours near the Fermi energy are concave and 
thermal population of concave isoen ergeti c contours can 
be neglected. Obviously, Av in Eq. (3.12) shows a criti- 
cal behavior: it is zero on the convex side and non-zero 
on the concave side of the transition. However, there are 



two other quantities which also show a critical behav- 
ior. As Figs, [sja) and[5jb) illustrate, even a concave FS 
does not necessarily have more than two self-intersection 
points: this happens only if the FS is shifted along cer- 
tain directions that lie close to high symmetry axes, i.e., 
q lies within some angular interval A^q, and the magni- 
tude of the shift is below certain threshold, i.e., q < q, 



Obviously, A6'q and gmax also depend on A in a critical 
manner. "^^ Approximating J cPq by lS.9qq^^^, we resolve 
the 5 functions and integrate over all energies to obtain 



5a,,, = ^ 2]^^ql^^k,.p„(qmax)|^ 
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(3.28) 



where the sum runs over all intersection points, the prime 
denotes a derivative with respect to the polar angle, and 
1 = 1/|1|. The task at hand now is to find the energy 
dependences of Av' and A6'q. 

This task is facilitated by the geometrical construction 
in Fig. |6f^b). Let 9* be the angle between the normal 
to the FS at a point, parameterized by the angle 9, and 
q. As one goes around the FS contour, 9* changes with 
9. Figure IgFc) shows the dependence of 9* on 9 for the 
FS in Fig. mh) in the convex (sp < Sc, dotted), crit- 
ical {sp = £c, dashed), and concave {ep > £c, solid) 



regimes. The dependence is monotonic for the convex FS 
and non-monotonic for the concave one. (This behavior 
is not specific to the particular FS considered here but 
is a general feature of any convex or concave contours). 
The oscillations are related to the rotational symmetry 
of the FS (six-fold in our case; Fig. ^c) shows only the 
domain 9 e [0,7r]). The non-monotonic parts are cen- 
tered around special ("invariant") points passed by the 
9* (9) curves for all types of contours. Near the invari- 
ant points, the non-monotonic part of the curve obeys a 
cubic equation 



9* = b9^ - a{A)9, 



(3.29) 



where a(A) oc A and 6 > is a constant. The energy 
dependences of the critical quantities can be obtained 
from this equation. 

To find A^q, we note that the equation Sk — £k-q = 
reduces to Vk • q = for small q. This implies that the 
solutions are those points on the FS where the normal to 
the FS is perpendicular to q [cf. Fig. letb)], i.e., 9*{9) = 



9q + 7r/2, where Sq is the angle defining the direction 
of q. From symmetry, if is a solution, so is 6* -f tt; 
therefore, one needs to consider only half the domain 
of 9. That only certain directions of q allow for more 
than one solution to this equation, may be appreciated by 
inspecting Fig. , which makes it obvious that multiple 
roots can only occur in the regions of non-monotonicity. 
The interval A9q where it happens is then proportional to 
the (vertical) width of these regions. Using Eq. (3.29), we 



find that A6lq oc A^* oc A^/^. Similarly, one can 



that Av* oc QmaxA and q„ 



oc 



A^/^ (a posteriori, this 



justifies the assumption of small q). 

Substituting these results into the expression for the 



max- conductivity, we find that A0q[Avj]^ oc A^/^, which 



means the prefactor of the term in the resistivity 
scales as A^/^. The T'^ term is always present, as dis- 
cussed before. Hence, the resistivity has the following 
form: 



p = Po + A 



-) 



9/2 



rp4 



(3.30) 



where po is the residual resistivity, 9{x) is the step func- 
tion, and A and B are material-dependent parameters 
(generically, A ^ B). A crossover between the T** and 
regimes occurs at T Sp^A/ep)^^^ <^ ep. 
Returning to the case of A < T, when both convex and 
concave contours are populated, it is easy to see that the 
^9/2 pi-efactor is replaced by T^/^, leading to a T^'^/^ 
term in p. This term, however, is subleading to the 



T one. Therefore, Eq. (3.30) describes the leading T- 



dependence of the resistivity in both situations (|A| ^ T 
and |A| < T) near the transition. Note that the expo- 



nents of 2, 4, and 9/2 in Eq. (3.30) are universal, i.e., 
they are the same for an arbitrary 2D Fermi surface with 
a non-quadratic energy spectrum near a convex-concave 
transition. 
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FIG. 6. (color online), (a) Isoenergetic contours for the 
surface states of the Bi2Te3 family of 3D topological insu- 
lators. The dashed line corresponds to the critical energy for 
the convex-concave transition, (b) For small q, points (black 
dots), where the normal to the FS is perpendicular to q, are 
the points of self-intersection, (c) 9* vs 9 [as defined in panel 
(b)]. Solid: ef > £c; dashed: ef ~ Ec, dotted: ef < £c;. (d) 
A zoom of the non-monotonic part of the graph in panel (c) 
Reproduced from Ref. 1141 courtesy of the APS. 



IV. HIGH-TEMPERATURE LIMIT 

So far, our analysis has been focused on the low- 
temperature limit, when the ee contribution to the resis- 
tivity is a correction to the ei one. From the experimen- 
tal point of view, however, it is important to understand 
whether the ee contribution may become larger than the 
ei one. It is the case for Umklapp scattering, whose con- 
tribution grows unabated up to the temperatures com- 
parable to the Fermi energy. The normal contribution, 
however, is different: it saturates in the limit when the 
ee relaxation time becomes shorter then the ei one. The 
effect of saturation was understood al ready in the earlier 
days of the electron transport theoryl^^SH very frequent 
ee collision establish a quasi-equilibrium state with the 
drift velocity fixed by ei scattering. The previous anal- 
ysis was, however, limited to the case when normal ee 
collisions affect the resis tivit y via the energy dependence 
of the e« relaxation time.'^'^In the next Section, we show 
that the saturation occurs even if the ei relaxation time 
does not depend on energy. 



can be viewed as a linear operator /ee acting on the non- 
equilibrium part of the distribution function /j^^^ = /k — 



/ee[/W](k)^^/ee(k,k')/^^', 



(4.1) 



The non-Hermitian matrix operator I^c can be repre- 
sented in terms of its left, <I>'^, and right, eigenstates 
as 



(4.2) 



where t*^ is the effective ee scattering time. The right 
and left states constitute an orthonormal basis: 



where V is the system volume in the D-dimcnsional space. 
A general solution of Eq. ( 2.19 ) can be expanded over the 
complete basis as 



(4.4) 



substituting this form into Eq. (2.19), we obtain an equa- 
tion for the coefficients cx: 



1 A 



(4.5) 



Only the zero mode (A = 0) contribution survives in the 
limit of 1/t*j, — oo, so the solution in the high-T regime 
is given by 



T-i-oo 



eria>^=0(k)^ E ^^=°(k')vk' • E<(4.6) 



It is not difficult to see that the right and left zero modes 
of are 



$,^=°(k)) = -ahn'^ ; (1>r"(k)| = Ch. (4.7) 



Indeed, 



(l>^="|/ec(xEk/cc(k,k')=0 



(4.8) 



due to momentum conservation; while 

/ec|$^=") CX E^cc(k,k')k'<, = (4.9) 



A. Saturation of the resistivity in a single-band 
metal 



We adopt the simplest model of point-like impurities 
with energy-independent scattering time, when the BE 
is given by Eq. (2.19). The ee collision integral in (2.19) 



because the collision integral, evaluated for equilibrium 
distribution functions, remains to be equal to zero if all 
energies are shifted as £k — >■ + u • k with u being an 
arbitrary k-independent vector. The zero modes form a 
_D-dimensional subspace labeled by the Cartesian indices 
i = 1,2, ... D [ summation over these indices is implied 
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modes (4.7 1 is 



in (4.4) and (4.6)]. The scalar product (4.3) of the zero- 



(*-="l<f-=°) 



C,C,y{EF){hk, 



where 



k 



(4.10) 



(4.11) 



The scalar product (4.10) is diagonal in the coordi- 



nate system associated with the principal axes of the 
quadratic form (kikj); normalization is ensured by choos- 
ing CiCi — [!/{ Ep){k^)]^^ . Using these properties, we 
reduce Eq. (|4.6|) to 



i k' 

nE^^(-<)W£^.- (4-12) 



-eTi 



Finally, we obtain the conductivity tensor in the high-T 
limit as 



(kf) 



(4.13) 



In the opposite limit of low temperatures, the standard 
expression reads 



2e'^v{EF)Ti{viVj) 



In contrast to the case of energy-dependent rj, when 
the low- and high-temperature limits of the conductivity 



differ in how Ti is averaged over the energy,'^'^ these lim- 
its in our case differ in how the conductivity is averaged 
over the FS. Naturally, the two limits coincide for the 
Galilean-invariant case. Notice that saturation holds for 
any dimensionality and shape of the FS, i.e., regardless of 
whether the temperature dependences of the resistivity 
starts with a or term at low temperatures, it will 
saturate at high temperatures. In reality, of course, other 
scattering mechanisms, such as electron-phonon scatter- 
ing, will mask the resistivity saturation. 

If the FS is not abnormally anisotropic, the low- and 
high-T limits are of the same order, which means that a 
true r^-scaling regime does not have room to develop. A 
mechanism in which such a regime is possible is consid- 
ered in the next Section. 



B. Two-band model: Scaling regime 

In this Section, we consider a simple model of a two- 
band metal with impurities. Since normal ee collisions 
affect the resistivity for any multiply-connected FS, we 
consider the simplest case of two bands with quadratic 
dispersions, s^^^ = /c^/2mi^2, and, in general, different 
impurity scattering times, tu and Ti2. We consider only 
the inter-band interaction (the intra-band one drops out 
in this case anyway) and neglect processes in which elec- 
trons are transferred from one band to another. The BE 
for this model in 2D can be solved exactly by generalizing 
the method of Appel and Overhauseil^see Appendix [a| 
with the result 



(4.14) p{T) 
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(4.15) 



where 



y/miTn2 
2Te% 



dqq 



W{q,uj) 



-co Jo Svr ^1 - (q/kFir^l - {q/2kF2f 



Lo^N (w) [N (lu) 



(4.16) 



and fc^"" 



min{fci?i, fci?2}- The result for the resistivity 
follows already from the equations of motiorP 



mi—— = -eE 
at 



1712- 



dvi 
'df 



-eE - 



mivi 

Til 
m2V2 

Ti2 



- ?7'^2(vi - V2) = 
-77ni(vi-V2) =0,(4.17) 



if the phenomenological "friction coefficient" 77 is ex- 
pressed via the microscopic scattering time as 77 = 
7r/e_FToc(T). 

An interesting case is when the masses are significantly 
different (as would be the case for a metal with partially 
occupied s and d bandi^, e.g., m2 ^ mi and conse- 



quently, Ti2 <C Til. At T — > 0, the two bands conduct in 
parallel, and the total resistivity is dominated by that of 
the lighter band 



P(0) 



1 



e^EF Til + Ti2 



e^EFTii ' 



(4.18) 



At T — > 00, the resistivity saturates at a value determined 
by the resistivity of the heavy band 
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TT 
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e^EF \ Til V"*2 
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e^£FTi2 



(4.19) 
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Therefore, the T ■ 
icantly 



and T = oo Umits now differ signif- 



p{oo) 
piO) 



Ti2 



7712 
7771 



> 1. 



The scahng regime, in which 



piT) 



TT 7772 



1 



e2 7771 epTcciT) 



(4.20) 



(4.21) 



occurs in a wide temperature interval Ti <^ T <^ T^, the 
boundaries of which are defined by 



1 



1 mi 
{Ti) Til 7772' 



1 



1 7771 
Ti2 m2 ' 



(4.22) 



This model can also be applied to the QPT, in which 
case it is natural to assume that critical fluctuations oc- 
cur only in the heavy band. Consequently, the effec- 
tive interaction is obtained from Eq. (2.6 1 by replacing 
vf — » vf2 and — > ?772/27r. 
(|4.16[) for (5 = and q < kf' 



Computing the integral 
we find the effective scat- 



tering rate in the NFL regime 



1 



167r3C(4/3) Tmlw^ 1 ^7^n2/3 



,{T) 8ir(2/3) 



7T7|a2 V 



y4/3 



(4.23) 

where Ci^) and r(a;) are the Riemann and F-functions, 
correspondingly. In this scenario, normal ee collisions do 
lead to a real scaling regime in the resistivity with an 
exponent given by "naive" power-counting argument. 



V. LIMITATIONS OF THE 
BOLTZMANN-EQUATION APPROACH 

The semiclassical BE does not capture two types of 
effects. The first type-quantum- results from quantum 
interference between ee and ei scattering; the second one 
-classical-from correlations in the electron flow patterns 
produced by different impurities. In this section, we dis- 
cuss the limits of validity of the semiclassical approach 
focusing on the 2D case. 



A. Quantum- interference effects 



scale that differentiates between low and high tempera- 
tures, i.e., between the diffusive and ballistic regimes, is 
Tdb = ^/t-i- For T <ti TdBi one is in the diffusive limit, 
characterized by a logarithmically divergent Altshuler- 
Aronov correction!^ with all coupling constants being 
of order one, \Sa\/aD ^ ln(l/TTi)/£i?Ti, where ao = 
e^EpTi/ir is the Drude conductivity. For T ^ Tdb: one 
is in the ballistic limit, where the correction scales lin- 
early with T: \5a\/aD ^ T/ep (Ref.lHl). Apart from the 
interaction correction, there is a also a weak-localization 
correction —SaY^n^/au ^ hi{T^/Ti)/£pTi, where is the 
phase-breaking time, a precise form of which depends 
on whether one is in the diffusive or ballistic limits: 
in the former, l/r^ ^ Tln{spTi)/epTi; in the latter, 
l/r^ ^ 1/tco- In the diffusive limit, the weak-localization 
correction is similar to the Altshuler-Aronov result, dif- 
fering only in the prefactor. In the ballistic limit, the 
weak localization correction is smaller than the inter- 
action correction by a factor of ln{T* /T)/TTi. There- 
fore, the correct order of magnitude for the quantum- 
interference correction is still given by the interaction 
correction both in the diffusive and ballistic limits. Com- 
paring the FL-contribution —Sapi^/ap) ^ T'^Ti/ep to 
the quantum corrections, we find that |(^(Tqc/(5(Tfl| 
\n{l /TT,)/T^Tf > I and \6aFi^/SaQc\ - Tr, > I in 
the diffusive and ballistic limits, correspondingly. There- 
fore, it is meaningful to consider the FL contribution and 
neglect quantum-interfence processes in the ballistic but 
not in the diffusive limit. The interplay of different mech- 
anisms is shown schematically in Fig. [7j 

In the integrable case, the term in the resistivity 
vanishes and the FL correction scales as \SaFi,\/aD ^ 
T'^T-Jsp. In this case, the FL correction dominates over 
the quantum one only at temperatures well above the 
diffusion-ballistic crossover: T » (eFTi)^/^TDB > Tub- 
It is worth noting that the discussion of the experimen- 
tal observations of quantum corrections in the ballistic 
regime has been so far limited to 2D electron gases in 
Si and GaAs heterostructures'*'*, with essentially circular 
FSs and almost parabolic spectra, where the FL contri- 
bution is expected to be very small. The FL contribu- 
tion, however, is expected to play a dominant role in 2D 
systems with highly anisotropic FSs and non-parabolic 
spectra, such as the surface state in the topological insu- 
lators of the Bi2Tc3 family, discussed in Sec. |IIIE 2| 



1. Fermi-hqmd regime 

Recall that the FL-like contribution to the resistivity, 
discussed in this paper, behaves as (or T**, if there 
is approximate integrability) in the low-temperature 
regime, defined by 1/tcc <C I/ti and saturates in the 
high-temperature regime, defined by I/t^c ^ l/'H- In 



non-integrable systems, I/tc, 



gT /ep, where g is 



the dimensionless coupling constant. In a generic FL, 
g ^ 1 and the crossover between the two limits occurs at 
T* ~ \/£p/Ti- For a good metal, epTi I so that 1/t,; <C 
T* <^ Ep. In case of quantum corrections (QC), the 



2. N on- Fermi-liquid regime 

In this section, we describe the interplay between the 
quantum-interference and direct ee contributions to the 
resistivity in the NFL-regime of a ferromagnetic quan- 
tum phase. For simplicity, we assume that the integra- 
bility is broken already in a single-band case by suffi- 
ciently strong concavity of the FS. In this case, one can 
simply calculate the transport time for scattering at crit- 
ical spin fluctuations, described by the propagator (2.6) 
with (5 = 0, and substitute the result into the Drude 
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FIG. 7. Different temperatue regimes for Fermi-liquid (FL) 
and quantum-interference (QC) corrections to the conductiv- 
ity. The shaded region on the temperature scale is the regime 
where the FL {T^) correction is dominant. 



formula.E2l This gives -(Jctnfl/c^d = T^/^ti/Tq^^, where 
To ~ (vp I a){kpaY (we remind that kpa 3> 1 is a con- 
trol parameter of the HMM model). The temperature 
above which the NFL contribution saturates is now given 
by T* - £^/*(fcFa)VTf/^ The main difference between 
the FL- and NFL-regimes is that quantum criticality 
changes space-time (or energy-momentum) scaling: while 
w oc g in a FL, uj (x near QCP, where Z is the dy- 
namical exponent (Z = 3 in the HMM model). There- 
fore, the temperature of the diffusive-ballistic crossover, 
determined by the condition qvpTi ~ 1, where q ~ 
{T/vpo^Y^"^ is a typical value of the momentum trans- 
fer in an ee collision, is replaced by Tdb ^ {a / v pTi)'^ / ti 
fRef. |43|) . In a clean system, where a/vpTi ^ 1, Tdb is 
significantly smaller than in a FL, where Tdb ^ ^/tii so 
that the ballistic regime continues down to much lower 
temperatures compared to the FL case, and we limit our 
analysis to this regime. (Lowering of Tdb in the vicin- 
ity of a ferromagnetic QCP becomes noticeable already 
in the FL regime.l^^l) Another consequence of quantum- 
critical scaling is that the quantum correction in the bal- 
listic regime near a QCP behaves as T^/^ as opposed to 
T: -SaQc/cTD - {T / epf'^ {akpY'"^ (Ref. IMl). Com- 
paring the NFL and QC contributions, we find that the 
NFL contribution dominates only at T ^ T, where 
f - {akp)^/Ti - TT:,Yi{spT;f{kpaf > Tdb- A relative 
weakness of of the NFL contribution is due to small-angle 
scattering at long-wavelength critical fluctuations. 



the ee interaction, l^c = vpTcc, is smaller than at least the 
average distance between impurities, , where Ni 

is the number density of impurities .ESln this regime, the 

ei mean free path k ^ ^ l^c is the largest scale 

of the problem, which implies that the FL contribution- 
even if allowed due to anisotropy of the FS-has already 
saturated of at a value comparable to the Drude resistvity 
(cf. Sec. IV I. Barring phonons, the viscous contribution 
is the only source of the T dependence in this regime. 

To estimate the magnitude of the viscous contribution, 
we consider, following Ref.|3Sl a flow of the electron liquid 
through a random array of spherical impurities. First, 
the impurity radius R is assumed to be much larger than 
lee, SO that a hydrodynamic description is applicable at 
all lengthscales. The force on one electron from all im- 
purities is just the Stokes force Fg ~ {Ni/N)imR, where 
/i is the dynamic viscosity, u is the flow velocity, and 
N is the electron number density. [An exact value of the 
numerical coefhcient in Fs depends on the boundary con- 
ditions for the velocity at the surface of the sphere^ but 
will not be needed here.] In steady state, Fg = eE, which 
yields u eEN /Ni^R, and thus the viscous contribution 
to the resistivity is given by 5pv ^ NifiR/e^N^. In a FL, 
5py ex. fjL ^ mNvplcc oc 1/T^; thus the viscous correction 
is of the insulating sign. On the other hand, the Drude 
resistivity resulting from scattering off the same impu- 
rities is pp) — m/e^NTi ^ mvpNiR^ /e^N . The viscous 
contribution is smaller than the Drude resistivity within 
the hydrodynamic regime: 5py/ pp) ^ lee/ R ^ 1- 

In 2D, the Stokes force from a disk-like impurity de- 
pends on R only logarithmically: Fg = Anpu/L, where 
L = \n{3.70p,/ RmNu) '.'^^ Rowever, l/n ~ vpNiR also 
contains R instead of R^, so that the ratio 5pv/ po is the 
same (up to a logarithm) as in 3D. 

The situation is somewhat different for small impuri- 
ties, because a force exerted by a small {R <C lee) sphere 
on a rarified gas depends not on the viscosity but on the 
gas-solid accommodation coefficients.*® In 2D, the sit- 
uation is further complicated vy the Stokes paradox.^S 
Hruska and Spivak**^ showed that the viscous correction 
for small impurities in 2D is given by 



5pv 
Pd 



le 



In 



(5.1) 



B. Viscous contribution to the resistivity 



The statement that ee interaction does not contribute 



to the resistivity of a Galilean-invariant FL (cf. Sec. II D 



) seems to contradict an intuitive notion that it is the vis- 
cosity of a liquid that defines its rate of flow. In a certain 
regime, indeed, the resisti vity d oes depend on the vis- 
cosity of the electron liquid.'^^'^ This effect is not taken 
into account by the standard BE which neglects not only 
quantum but also classical correlations between scatter- 
ing events. The "viscous" contribution occurs at high 
enough temperatures, when the mean free path due to 



where a is the impurity scattering length and a <C Zee 

1 /2 

l/A^j . Because of a large logarithmic factor, 5py can, 
in principle, be comparable io po- 



VI. 



CONCLUSIONS 



The main purpose of this paper was to analyze the ef- 
fect of ee interactions on the resistivity in the situation 
when Umklapp scattering of electrons can be neglected. 
Such a situation arises, e.g., in low-carrier density mate- 
rials, as well as in metals near a g = QPT, where the 
effective interaction is of a long range. In such cases, the 
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conventional dependence (or its analog in the 

NFL region near a QPT) of the resistivity on tempera- 
ture is not guaranteed. Whether it is present depends on 
1) dimensionality, 2) shape, and 3) topology of the FS. 
If the FS is quadratic or isotropic, there is no con- 
tribution to the resistivity. However, anisotropy is not 
sufficient to guarantee the dependence. In the case 
of a convex and simply connected FS in 2D, there is no 
dependence either. In such cases, the leading tem- 
perature dependence on resistivity due to ee interactions 
is in the FL region and t(^+8)/3 in the NFL region. 
Also, if the FS changes its shape from convex to con- 
cave as a function of the filling fraction, the resistivity 
follows a universal scaling form near the convex-concave 
transition. In all other cases, the r2 (or T^) behav- 
ior is allowed, albeit only as a correction to the Drude 
resistivity. However, a true scaling regime (when the ee 
contribution is larger then the ei one) is possible for a 
quantum-critical two-band metal with substantially dif- 
ferent band masses. Since a quantum-critical behavior 
is observed typically in multi-band metals with partially 
occupied d bands, we conjecture that the 5/3 scaling of 
the resistivity observed in 3D ferromagnetd^ and sub- 
quadratic scaling in a quasi-2D metamagnet Sr3Ru207E2l 
is due to the interaction between light and heavy carriers 
in these materials. 
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Appendix A: Two-band model in 2D 



In this Appendix we derive Eq. (4.15). The two coupled BEs (1 and 2 refers to the bands 1 and 2) read 

/i(ki) 



evi • Vki/i 
-ev2 • Vki/2 



[/l,/2] 



T\2 



(Al) 



where PJ^ is the electron-electron collision integral for scattering between two electrons from the ith and jth band, 
and 2 are the equilibrium distributions. Since in our model each of the bands is Galilean-invariant on its own, the 
intra-band ee interaction cannot affect the resistivity, and the corresponding parts of the collision integrals are not 
written down. Linearizing (a = 1, 2) in the same way as for the single-band case 



= + Ua (1 - na)go 



we obtain 



evi-En'j 
-ev2-En2 



-lie [91:92] 
-Ice [91,92] 



Tn[gi 

Til 
T 71292 
T\2 



(A2) 



(A3) 



where 



rl2 



f ^W^(k-k',ek,-£kOni(ki)n2(k2)(l-ni(k;))(l-n2(ki)) 
X [51 (ki) + 92 (k2) - 91 (k'l) - 92 (K)] ^ (^ki + £1.2 - £k[ - S (ki+ka-k^ - k'2) 



(A4) 



and I^^ differs from I^^ in that the first integral goes over ki instead of k2 . We seek for a solution in the following 
form 



ga = eca (vq • E) /T, 



(A5) 
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where Cq are the constants that are to be determined. Multiplying the first and second BEs by Vi • E and V2 • E and 
integrating over ki and k2, correspondingly, we obtain for the left-hand sides 



i^c ■ E) 



{2ny ^ " ' " 47r 
Similarly, the ei collision integrals, integrated over the corresponding momenta, reduce to 



(27r)^ Tia 



Va-E = 



An 



(A6) 



(A7) 



For a FL, the scattering probability may be assumed to depend only on the momentum transfer but not energy 
transfer. In this case, the integrals of ee collision integrals multiplied over • E reduce to 

''^\v, . E/- = eE^R- 



where 



(2»)= 



VF1VF2 



mi 



_C2_ 



2k" 



dqq 



W{q) 



\2kFi J 



(A8) 



(A9) 



with fc™'" = min{fci?i, fci^2}, and V2 • E/^^ differs from ( A8) by a factor of mi/m2- [Integration over energies was 
performed with the help of Eq. (|3.7[).] Solving the system of linear equations ci and C2 



mi 

m2 
47r 



Vf2 = 



ci mi 
Til 47r 
C2 m2 
Ti2 47r 
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m2 



R 



VfiVF2 
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VF1VF2 
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m2 
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(AlO) 



we find 



Cl 



C2 



1/Ti2 + (1/tcc) {mi/m2 + 1) 
1 I ]_ [ JL I L "12 ) 

T-iiri2 Tco \Tii m2 ri2 mi ^ 
1/rii + (1/Tee) (m2/mi + 1) 



1 1 mi 

Tco \Til m2 



1 m2 
Ti2 mi 



where the effective ee scattering time was introduced as 

1 

Tr 



AttR- 



^m\m2 Vp^v 



(All) 



(A12) 



F2 



Using that the Fermi energy is same for both bands, i.e., that vfi/vf2 = Y^m2/mi, we cast Eq. ( A12 1 into a different 
form 



1 



= — T- 
3 



2k" 



dqq 



W{q,0) 



^l-{ql2kFiY^l-{q/2kF2Y 



(A13) 



Once ci9 are found, one readily finds the electric current and arrives at the expression for the resistivity quoted in 
Eq. (4.15). An explicit expression for I/tcc was derived in Ref. i50i for a special case of the " overscreened" Coulomb 
potential. Our result for this case coincides with that in Ref. 1501 up to a numerical coefficient. 

In general, the ee scattering probability depends not only on q but also on u. In this case, the integrals over 
energies cannot be performed in a general form, and the expression for I/too can only be reduced to the form quoted 
in Eq. (4.16). It should be stressed that the electron masses occurring in all equations of this section should be 
understood as hare rather than renormalized masses. This follows from the derivation of the BE in the Keldysh 
technique using the Migdal-Eliashberg approximation. 
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